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Abstract. We consider the shape evolution of two-dimensional islands on a 
crystal surface in the regime where mass transport is exclusively along the is- 
land edge. A directed mass current due to surface electromigration causes the 
island to migrate in the direction of the force. Stationary shapes in the pres- 
ence of an anisotropic edge mobility can be computed analytically when the 
capillary effects of the line tension of the island edge are neglected, and con- 
ditions for the existence of non-singular stationary shapes can be formulated. 
In particular, we analyse the dependence of the direction of island migration 
on the relative orientation of the electric field to the crystal anisotropy, and 
we show that no stationary shapes exist when the number of symmetry axes 
is odd. The full problem including line tension is solved by time-dependent 
numerical integration of the sharp-interface model. In addition to stationary 
shapes and shape instability leading to island breakup, we also find a regime 
where the shape displays periodic oscillations. 
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1. Introduction 

The interest of researchers in the phenomenon of electromigration has in the past 
mainly been guided by the technological implications of this effect, namely the 
degeneration of integrated circuits ^|2H2| . A large body of experimental work has 
been devoted to the study of electromigration- induced failure of metallic conductor 
lines under accelerated testing conditions OEJODEI- At least for simple geometries 
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(e.g., a single void in a single crystal grain), such experiments have been success- 
fully modeled using a macroscopic continuum theory of shape evolution [SJ El (see 
Sect|2J|. On the other hand, theorists have conducted a lively (and still ongoing) 
debate on the microscopic driving force of electromigration . 

What appears to be largely missing in the field is experimental and theoret- 
ical work on the mesocopic level, addressing the evolution of simple atomic-scale 
structures such as individual steps and single-layer islands. This could help to 
bridge the gap between the elementary moves of single atoms and the failure be- 
havior of polycrystalline conductor lines with a complex internal geometry. In the 
present paper we address the electromigration-induced shape evolution of single, 
atomic-height islands on a crystal surface. We work in a continuum setting first 
introduced, in this context, by Pierre-Louis and Einstein |l()j . Among the different 
kinetic regimes considered in their work, we focus here on the particularly simple 
case where mass transport is exclusively along the island edge, so that the area of 
the island is strictly conserved and the evolution law for the island boundary is lo- 
cal. Our new contribution consists in including crystal anisotropy, which is clearly 
necessary to make contact with island electromigration on real crystal surfaces 
(such as silicon |11M12| 1 and in lattice Monte Carlo simulations [11)11131 HI) . 

The continuum model is introduced in the next section, and previous work 
on the problem is described. In Section |3] we show how stationary shapes can 
be explicitly computed when the smoothcning effects of the edge free energy is 
neglected, and numerical results for the full, time-dependent problem are presented 
in Section 0] 

2. Continuum Model of Shape Evolution 

The physical system consists of an island of monoatomic height on an otherwise 
flat crystal surface in interaction with the flow of the electrons. The island should 
be large enough to allow for a coarse grained description where the individual 
particles of the assembly are blurred into a continuous entity. To put this into a 
mathematical setting we represent the spatial conformation by a closed curve in the 
xy-plane parametrised by the arclength s. All geometrical and physical quantities 
can now be expressed as functions of s. We describe the local orientation of the 
island edge by the angle 9 between the normal and the y-axis (counted positive 
in the clockwise direction). The derivative of 9 with respect to s is then the local 
curvature k. The force F acting on an atom at the island edge is composed of two 
contributions: 

F=-^n) + F el . (2.1) 

Here 7K is the chemical potential of the island edge 1 , with 7 denoting the edge 
stiffness, which is related to the step free energy 7(0) by 7 = r y + d 2 j/d9 2 [To] , The 
gradient of the chemical potential on the right hand side of (|2.1|l accounts for the 



1 Throughout we work in units where the atomic area and the lattice spacing are set to unity. 
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effect of capillarity which tends to minimize the free energy by driving material 
away from regions of high curvature. The second term is the force that actually 
causes electromigration. It is conventionally written in the form F e i = eZ*E t , 
where e is the elementary charge, E t is the component of the local electric field 
which acts tangentially to the step edge, and Z* is the effective valence of the 
step atom. On this level of modeling, all microscopic aspects of electromigration 
are lumped into the quantity Z*. To give an example, the effective valence for a 
copper atom moving along a close-packed step on the Cu(100) surface has been 
calculated to be Z* w —31 |14j . The form of the tangential electric field to be used 
in this work will be specified later. 

The motion of the atoms along the edge can now be described by a mass 
current j that is proportional to the acting force, with the factor of proportionality 
defining the orientation-dependent step edge mobility o~(Q), 

3 = °(P)F. (2.2) 

If we assume that atoms can neither detach from the island nor attach to it, 
the equation of shape evolution is simply the continuity equation that relates the 
divergence of the mass current to the normal velocity v n of the island edge, as 

Vn = ~h = " eZ * Et] - (2 - 3) 

By dimensional analysis, the comparison of the two terms inside the square brack- 
ets on the right hand side of 1)2.3)1 defines the characteristic length scale |l()lll(jllT7j 



'° - (2A) 

where Eq is a typical value of the electric field strength. For islands small compared 
to I e the capillary forces dominate the evolution, while for islands large compared 
to I e the evolution is dominated by the electromigration force. Thus Ie sets the 
length scale on which electromigration-induced shape instabilities can be expected. 

As formulated so far, the model is identical to the one that has been used 
extensively to model the evolution of quasi-two-dimensional (cylindrical) voids in 
current-carrying metallic thin films joT m IT51 ITT?1 r^T)l PT1 [^1 PTT1 f^"TT IUjI f?ol f^Hl [T71 
128) . In that context it is essential to take into account the effect of current crowding, 
which refers to the fact that current lines cannot pass through the insulating void, 
and hence the electric field is enhanced in the remaining parts of the conductor. 
As it is necessary to monitor the changes in the electric current configuration 
that occur in response to the shape changes of the void, the dynamical evolution 
becomes manifestly nonlocal. In the physical situation of interest in the present 
paper, this complication does not arise, since a single-layer island is not expected to 
significantly modify the electric current configuration in the underlying substrate 
film. We may therefore assume a constant electric field of strength Eq directed at 
an angle cj> with respect to the x-axis, and set E t = Eq cos(9 + <f>) in l|2.3|l . 
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A model that interpolates between the case of a completely insulating void 
and the case of a constant electric field considered here is obtained by assign- 
ing a conductivity S' to the void which differs from the conductivity E of the 
bulk material The two limiting cases then correspond to E'/E = and 

E'/E = 1. The local model with E'/E = 1 has also been used to describe the 
electromigration-driven evolution of dislocation loops |291 130j . 

Our primary aim in this paper will be to determine the stationary shapes 
of the drifting island. The condition for an island to move at constant velocity V 
along the cc-axis without changing its shape is 

v n = Vsm{9). (2.5) 

It has been shown long ago that in the absence of anisotropy (7 and a independent 
of 6) a circle is a stationary solution of l|2.3|l . for general conductivity ratio E'/E 
[T%| . In the following we try to find corresponding results for the anisotropic case. 



3. Analysis in the Absence of Capillarity 

Combining the equation of motion (|2.3|l and the stationarity condition (|2.5|l gives 
d d 

Q- a °W)\Q- s (7«) - ^0 cos(0 + 0)] = Vsin(0) (3.1) 

with Fq = eZ* Eq. This is a nonlinear ordinary differential equation of fourth order 
for the shape, or alternatively of second order, when viewed as an equation for the 
curvature as a function of 9 (using = tj^)- In this section the capillarity term 
proportional to the stiffness 7 will be neglected. While it is true that the effects 
of capillarity become small for large islands (island radius 3> Ie), it can hardly 
be justified to drop this term altogether, because it contains the highest order 
derivative in the equation. In fact islands that are large compared to Ie break 
apart, rather than reaching a steadily drifting shape JOj (see Sect^J). 

Nevertheless the analysis of the capillarity-free case 7 = is rather instruc- 
tive. It is accessible to an analytical treatment, because the unknown curvature 
enters (13.111 only through the derivative Jj which can be eliminated by the follow- 
ing substitution: From the two geometrical relations 

{dsf = (dxf + {dy? «» (g) 2 + (g) 2 = 1 (3.2) 

and 

dv ^ 

tanW = -£ = --t (3 - 3) 

ds 

we conclude that || = — sin(#). Replacing sin(0) in l|3.1|) using this relation we 
find that the derivative of Vy — j is zero and therefore 29 Vy = j + const. 
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a{0) = UUL V ol^^T (3-6) 



The constant of integration can be set to zero, since it generates only a vertical 
displacement. Hence the y-coordinate of the island shape as a function of 9 is 

y(e) = ^-a(8)cos(0 + cj)) (3.4) 

and using = — tan(#) once more, we obtain 

x{9) = - j\ot{8')^d8'. (3.5) 

In order to derive stationary island shapes from l|3.4|l and l|3.5|l one has to 
specify the mobility u{9). To mimic the anisotropy of a crystal surface with a 
certain discrete rotational symmetry, we choose |16] 

g [l + 5'cos 2 (£>((9 + a))] 
l + 5cos 2 (£»a) 

where a is the angle between the symmetry axes and the coordinate system 2 , 2D 
is the number of symmetry axes, and S > is a measure for the strength of the 
anisotropy; see FigOfor an illustration of i|3.6[l . The denominator ensures that the 
mobility at 9 = is crn. 

3.1. Existence of stationary shapes 

When calculating the shape by means of (|3.4() and l|3.5|l with a given choice of 
D,S,a and </), one encounters two types of irregularities: (i) The curve is generally 
not closed, and (ii) it can contain self-intersections. To get a closed shape one has 
to demand that x and y are 27r-periodic functions of 9. Because y{9) is always 
27r-periodic it suffices to have 

/ dx{6)= d9 cot(0)[-£cos(fl + 0) - a sin(0 + </>)]= 0. (3.7) 
Jo Ja d9 

For integer values of D (even rotational symmetry) the mobility is 7r-periodic and 
hence the integrand in H3.7(l is odd under shifts 9 — > 8 + n. As a consequence (13.7(1 
is always fulfilled. For half-integer D (odd rotational symmetry) the evaluation of 
the integral in 1)3. 7fl yields 2D cos(0) cos(2Z)a) — sin(0) sin(2Da!) times a nonzero 
factor, so that l|3.7|l requires that 

tan(2Da) tan(0) = 2D. (3.8) 

The generation of self-intersections is illustrated in Fig^ Because the self- 
intersections are always accompanied by singularities (sign changes) of the curva- 
ture, they can be avoided by requiring that 

1 1 dy 



k sin(0) 88 



> (3.9) 



2 The reason why we do not choose a symmetry axis to coincide with one of the coordinate axes 
is that the stationarity condition 12.51 becomes particularly simple when the direction of motion 
coincides with the ^-direction. 
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S = 0.35 S = 2 S = 10 



Figure 1. Formation of self-intersections in the capillarity-free 
case. The figure shows stationary shapes for D = 3 and a — cp = 0, 
where the critical anisotropy strength is S c w 0.35. 

everywhere. Thus the smoothness condition is related to the existence of extrema 
of y{9). If || is nonzero at or n, the integral in (|3.5|) diverges because it contains 
the singularity of the cotangent function; the shape is then not bounded. If ^| is 
zero somewhere else (|3.9f) is violated. If D is half-integer, the condition §f (0) = 
§|-("7r) = requires a = <j> = 0, which violates (|3.8|) . In other words: For half-integer 
values of D the shape is either not closed or not bounded. While this result has 
been derived using the specific functional form (|3.6|) for the mobility, we expect 
it to be true in general that in the capillarity-free case, no stationary shapes exist 
when the number of symmetry axes is odd. 

In the following we may therefore restrict ourselves to anisotropies of even 
symmetry. In this case §f (0) = || (tt) because of the 7r-periodicity of y and we 
end up with two conditions, 

dy , , SDsm(2Da) , , . 

of) 1 + 6 cos z (Da) 

For a given value of the angle a — <j> between the symmetry axes and the field, 
(|3.10|) determines the direction of motion of the island, and the inequality l|3.11|) 
determines the values of S for which a smooth shape exists. For S = there is 
always such a shape for every choice of a — 4>, namely the circle which moves in 
the direction of the applied field. Increasing the anisotropy strength S the circle 
elongates along the field direction and finally at some critical value S c sharp corners 
develop as precursors of self-intersections (see FigQ. For S > S c no stationary 
shapes exist. Figure |3 shows the array of possible shapes for the case of sixfold 
symmetry (D = 3). Note that in the absence of capillarity the characteristic scale 
(|2.4(l disappears from the problem, and hence the shapes are independent of island 
size (the island size determines however the migration speed, see Sect. 

It is straightforward to derive the critical anisotropy strength S c for the 
special case a = n/2D, = 0, where the direction of island motion and the field 
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Figure 2. Stationary island shapes without capillarity, for six- 
fold crystal anisotropy and different values of the angle a — <j) 
between the symmetry axis of the anisotropy and the field, and 
the anisotropy strength S. In this figure the field direction is hor- 
izontal, and the arrows indicate the direction of island migration. 
In the empty part of the (S, a — 0)-plane no stationary shapes 
exist. 



direction coincide with the direction of minimal mobility. The island then elongates 
perpendicular to the field and the self-intersections are forced by symmetry to 
appear at 9 = and it. The condition of vanishing curvature radius reads simply 
d 2 y/d8 2 \ e=0 — 0, which implies [a — d 2 a/d9 2 ]\g =0 = 0. For our choice (|3.bjl of the 
mobility function this yields 

Sr= „J : for = 0. (3.12) 



2D 2 



1 



This is the minimal value of S c ; the maximum range of smooth shapes appears 
at a = <p = 0, where the island moves along the direction of maximal mobility 
(FigEJ. 

3.2. Direction of island migration 

An important consequence of anisotropy, which remains true also when capillary 
forces are turned on, is that the direction of island motion does not generally 
coincide with the direction of the electromigration force when a — <p 7^ 0. This 
effect was previously observed for the nonlocal model ^5] . The relationship i|3.10|) 
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shows that a and <f> have opposite signs, which implies that the direction of island 
motion lies between the field direction and the symmetry axis. 
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Figure 3. Relative angle 4> between the field direction and the 
direction of island motion as a function of the angle a — <j> between 
the field and the axis of maximal mobility, as computed from 
i|3.10[l with D = 3. The vertical line marks the minimal mobility 
direction a = ir/6. 

Figurc|31shows the angle between the direction of island motion and the field, 
cj>, as a function of a— <f>, which is the physical control parameter determined by the 
experimental setup. The island moves in the direction of the field, <fi = 0, both when 
the field is along the maximal mobility direction (a — <j) — 0) and along the minimal 
mobility direction (a — (f> = n/2D). With increasing anisotropy strength the graph 
becomes strongly skewed towards the minimal mobility direction until, at a second 
critical anisotropy strength given by S c — 1/2D 2 , the function becomes multi- 
valued. Beyond this point <fi approaches a nonzero limit as a — <fi — > tt/2D, despite 
the fact that cj> = exactly at a — <\> — tt/2D. The physical consequence is that 
the direction of island motion, as well as the island shape, change discontinuously 
as the field direction is moved across the direction of minimal mobility. Since S c 
is (slightly) smaller than the minimal value of S c given by (|3.12|l , there is a small 
range of anisotropy strengths where smooth shapes exist for all angles but the 
shape changes discontinuously at a — <fi — n/2D. 

3.3. Migration speed 

Since both coordinates of the parametrization l|3.4|l and l|3.5H are multiplied by 
1/V it is clear that the velocity is inversely proportional the extension of the 
island; for a circular island of radius R it is well known that V ~ 1/R JSj- I n 
other words, for every allowed choice of a — <j> and S, l|3.4|) and l|3.5|l generate a 
whole family of similar shapes moving with different velocities. To determine the 



Electromigration of Islands with Crystal Anisotropy 



dependence of V on the island size we calculate the area A of the island as follows: 
A = jj y{6) dx{9) = -\fj d9 cot(0)^y(9) 2 



d9 cot(0) — cr(6») 2 cos 2 (6»- 
dO 



_1 (Fq 

Inserting the form (|3.6|) for the mobility, it follows that the migration speed is of 
the form 

Fqctq ^aS 2 + bS + TT 

l + Scos 2 (Da)' [ ' 

where the coefficients a and b are integrals of combinations of trigonometric func- 
tions which depend on a, cf> and D. For the special case a = <j> = and D = 3 
illustrated in the bottom row of Fig|2 we find a — 397r/8 and b — 7ir, which 
implies that the speed increases by about a factor of 1.5 when the anisotropy is 
increased from S = to the critical value S c « 0.35. The speed increase is related 
to the elongation of the island shape, which brings the orientation of the edge 
closer to the maximum mobility direction 9 = 0. When the island moves along 
the direction of minimal mobility (a = tt/2D, = 0) the coefficients in Q3.13H are 
negative, a = —9ir/8 and b = — 57r, which implies that the island slows down with 
increasing anisotropy. 

4 . 
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4t0 
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Figure 4. Stationary shapes without anisotropy. The dimension- 
less initial radii of the different shapes are po = 3.3, 4, 5 and 6, 
and the computation was done with 100 discretization nodes. The 
dotted line shows the circular shape which constitutes the stable 
stationary solution for po < p c ~ 3.26. 



4. Numerical Integration of the Time-Dependent Problem 

We now integrate the full dynamical problem (|2.3|l by means of the sharp inter- 
face method described in |16l l3~T] . The evolution equations are discretized and 
integrated using a variable-step variable-order predictor-corrector method [22] • To 
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avoid frequent remeshing an additional tangential velocity was introduced which 
keeps the nodes equidistant and which has to be computed implicitly. 

Again we are mainly interested in the long time behavior of the system. As 
initial condition we take a circular island, usually imposed with a slight deforma- 
tion to trigger a certain mode of motion 3 . We restrict ourselves to the case of a 
mobility with six-fold symmetry (D = 3) with the field aligned parallel to a sym- 
metry axis. The stiffness 7 is assumed to be isotropic 4 . The tunable parameters 
are then the initial radius Rq and the strength S of the anisotropy. All lengths 
are expressed in units of the characteristic scale / e , and the dimensionless initial 
radius is pg = Rq/Ie- 

10.0 1 . 1 1 1 1 1 . 1 . 1 

5 ■ - y^Xy^X^^^S ) 



0.0- 



-5.0- 

10.0 - 

10.0 0.0 10.0 20.0 30.0 40.0 

Figure 5. Breakup of an island in the isotropic case. The figure 
shows a time sequence of shapes, with a shading that darkens 
with increasing time. The initial shape is an ellipse elongated in 
the direction of island motion, with an area corresponding to a 
circle of dimensionless radius po — 7. The computation used 150 
discretization nodes. The breaking of the up-down symmetry is 
due to numerical noise. 

To link our investigation to previous work we start with the isotropic case 
5 = 0. Here we find circular stationary solutions for small radii. At the crit- 
ical dimensionless radius p c — 3.26 the circular solution loses its stability [T§| 
and a bifurcation to non-circular stationary shapes occurs. In |30| this bifurcation 
was investigated by numerically solving the stationarity condition (|3.1fl . and two 
branches of post-bifurcation shapes were found. In our time-dependent integration 
we only observe one branch, referred to as 'mode I' in |30j . To linear order in 
the distance from the bifurcation, the second 'mode II' branch is the symmetrical 

3 Unless stated otherwise, for the examples described here the asymptotic mode of evolution is 
independent of the precise initial condition. 

4 This is physically motivated by the fact that the anisotropy of the mobility <r, which is a 
thermally activated quantity, is expected to be much more pronounced than that of 7. 
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counterpart with respect to the circle, in the sense that Ki(0) + Ku(0) — const.; 
mode I is elongated in the direction of the field, mode II perpendicular to it. The 
mode II shapes which move at a smaller velocity seem not to be selected by the 
dynamics. 

Some noncircular stationary shapes for S = are shown in Fig0] We see 
that for po > 3.5 the shapes become non-convex. These shapes were not found 
in |30| . because the numerical calculation was stopped at the point where 1/re 
approaches zero. Increasing the dimensionless island size further above po — 6.2 
causes the island to break. It first evolves a long finger, which corresponds to the 
'slit'-solution described in 17 , and which subsequently detaches to form a new 
island (Fig. |5j). Similar breakups occur in the presence of anisotropy, so that the 
range of po is generally restricted towards large values. 









) 




Figure 6. Left panel: Stationary shapes with anisotropic mobil- 
ity of sixfold symmetry (D — 3). The initial shape was a circle 
of dimensionless radius po = 2.5, and the anisotropy strength is 
S = 1, 5, 10 and 100. The computation used 50 discretization 
nodes. Right panel: Polar plot of the mobility corresponding to 
the shapes in the left panel. 

For the anisotropic case there are no circular stationary shapes anymore. 
The deviation from the circle increases with increasing size po and strength of 
anisotropy S. For strong anisotropy the islands evolve facet-like features, i.e. re- 
gions of the edge where the radius of curvature is very large. Figure El shows some 
examples. As was already discussed in Sect. 13.21 the direction of motion generally 
differs from the field direction, when the latter is not aligned with the symmetry 
axes of the anisotropy. Under certain conditions the symmetry of the direction 
of motion with respect to the field and anisotropy direction appears to be spon- 
taneously broken, i.e. we have observed stationary island migration off the field 
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direction even when the latter coincides with a symmetry axis of the anisotropy 
|33| . This behavior may be related to the discontinuity of the direction of mo- 
tion when the field direction crosses the direction of minimal mobility, which was 
described in Sect V3. 21 for the capillarity-free case. 

With increasing S the range of po for which stationary solutions exist shrinks 
toward smaller values and is replaced by a domain where a qualitatively new 
behavior appears. Here the time evolution neither runs into a steady state nor 
into a breakup event but instead becomes oscillatory. A typical sequence is shown 
in Figd First the island elongates along the field direction, then one side of it 
bends so that the conformation becomes nonconvex. This leads to an acceleration 
of the mass transport which then shrinks the shape back to its original state. 
The oscillatory shape evolution in Fig[7]is moreover seen to break the symmetry 
with respect to the field and anisotropy direction (recall that the field is directed 
along the x-axis which corresponds to a direction of maximal mobility, see the 
right panel in FigEJ ■ The solution shown in Fig[7| coexists with one obtained by 
reflection at the x-axis, in which the island travels to the right and down. Which of 
the two solutions is realized depends on the initial conditions. For other parameter 
values we have also observed solutions that switch between periods of upward and 
downward motion [2Sj. 

The conditions for the appearance of stationary and oscillatory shapes largely 
remain to be clarified. On a qualitative level, we have observed that the "facets" 
that dominate the stationary shapes for large S tend to be close to orientations 
of maximal stability in the sense of linear stability analysis [31]. In contrast, the 
oscillatory shapes are typically bounded by one stable and one unstable orientation, 
and the shape oscillation can be viewed as the propagation of a wave along the 
unstable side of the island (see Fig[7J). 

15.0 I . 1 ■ 1 . 1 . 1 ■ [ 



10.0- 
5.0 - 




0.0' 1 1 1 1 1 1 1 1 1 1 

10.0 20.0 30.0 40.0 50.0 60.0 

Figure 7. Oscillatory island motion. The dimensionless initial 
radius was po — 4 and the anisotropy strength S = 1. Note that 
the upper boundary is stable, while the lower boundary supports 
a traveling wave. The computation used 100 discretization nodes. 

For the reasons described in Sect [3] a quantitative comparison between the 
shapes computed numerically for the full problem and the analytically derived 
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capillarity-free shapes is not meaningful. On a qualitative level, the most apparent 
difference between the shapes in Fig[2] and Figs l4l6l is that in the presence of 
capillarity the direction of motion can clearly be discerned from the shape, while 
the capillarity- free shapes in Fig[2] are front/back symmetric. This can be seen 
directly from l|3.1fl . If 7 = the equation reduces to (again using = 

FK^-a{9) cos{9 + <j)) = V sm(9) (4.1) 
dO 

which possesses two symmetries: First, reversing the direction of the field and the 
velocity simultaneously leaves the equation invariant. The second symmetry is the 
7r-periodicity of the equation 5 which results in periodicity of the solution. Both 
symmetries are destroyed for nonzero stiffness. 



5. Conclusions 

It is interesting to compare the local (E'/E = 1) and nonlocal (E'/E < 1) boundary 
evolution problems in the light of our work 6 . For the isotropic case, it is known that 
the circular solution becomes linearly unstable for sufficiently large radius when 
E'/E > |19| . while it retains its stability (but becomes nonlinearly unstable under 
a sufficiently large perturbation I221EZ]) when E'/E = |191 \20\ . Furthermore it 
has been shown by complex variable techniques that noncircular stationary shapes 
cannot exist for E'/E = [221, hence a void is forced to break up when it becomes 
unstable (221 - Here we have shown that breakup occurs for sufficiently large islands 
also in the local case, but in addition there is an intermediate regime, between the 
linear instability of the circular solution and breakup, where noncircular stationary 
solutions are stable 7 . 

Perhaps our most remarkable finding is that crystal anisotropy may induce 
complex, oscillatory evolution modes in the local problem. An example of oscil- 
latory void evolution was previously reported by Gungor and Maroudas in the 
nonlocal case Future work should address the robustness and generality of 
the phenomenon. For this purpose it should be useful to employ numerical meth- 
ods that are complementary to the sharp interface algorithm adopted here, such 
as finite element [21], phase-field |23ll26ll25| . level set [21] and kinetic Monte Carlo 
[11)11131 IT4] approaches, and to extend the study of the additional kinetic regimes 
(other than the case of edge diffusion considered here) in Q2| to include crystal 
anisotropy. 



5 The mobility is 7r-periodic only for integer D, which is also the reason why we always get closed 
forms in this case (see Sect 1^1. 

6 For a similar comparison in the substrate geometry see |35|. 

7 Note that in I3UI the stability and dynamical relevance of the noncircular stationary solutions 
was not addressed. 
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